Determination of Landau Fermi-liquid parameters of strongly interacting fermions by 

means of a nonlinear scaling transformation 
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A nonlinear transformation approach is formulated for the correlated fermions' thermodynamics 
through a medium-scaling effective action. An auxiliary implicit variable-effective chemical potential 
is introduced to characterize the non-Gaussian fluctuations physics. By incorporating the nonlocal 
correlation effects, the achieved grand partition function is made of coupled highly nonlinear para- 
metric equations. Analytically, the low temperature expansions for the strongly interacting unitary 
Fermi gas are performed for the adiabatic compressibility-sound speed and specific heat with the 
Sommerfeld lemma. The expressions for the Landau Fermi-liquid parameters Fq and F{ of the 
strongly interacting fermion system are obtained. As a universal constant, the effective fermion 
mass ratio is m* /m = ^ at unitarity. 

Keywords: Rigorous results in statistical mechanics, Bose Einstein condensation (Theory), Series 
expansions 
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The Landau theory of Fermi liquid plays the key role 
in understanding the novel interacting fermions quantum 
many-body physics [lj. Either in the non-relativistic or 
relativistic frameworks, the central task is how to cal- 
culate the Landau Fermi-liquid parameters [2, y, |j, Q, 
which contributes to addressing the bulk or kinetic trans- 
port properties in terms of the quasi-particle viewpoint. 
The conventional techniques for calculating the effective 
fermion mass and/or ground state energy will involve the 
various coupled integral equations for the multi-points 
correlation Green functions exhibited by the Galitskii, 
Bcthe-Sepeter or Bcthe-Goldstone expansions [fj|. 

In recent years, considerable efforts for understand- 
ing the physics from Bardeen-Cooper-Schrieffer to Bose- 
Einstein condensation(BCS-BEC crossover) with ultra- 
cold atomic Fermi gases have been made. The two-body 
interaction strength can be tuned with the Feshbach reso- 
nance. At unitarity, the divergent scattering length with 
a zero energy bound state can manifest the universal 
properties[zfll©,[l(J. 

The fundamental issue on the unitary fermions ther- 
modynamics is on the symmetric fermions ground state 
energy [3, Q- Another exciting one concerns the asym- 
metric system thermodynamics. To describe or interpret 
the phase separation properties achieved by the experi- 
mentally trapped systems, the Landau effective fermion 
mass becomes the key dynamical parameter QJ], EH H3J . 
Like in the ground state energy, the infinite scattering 
length will certainly drop out in the effective mass expres- 
sion, which excites considerable attempts [HI. H2. [HI Ha] . 
The existing values are in the regime m* jra ~ 1.04 — 2.5. 
What is the exact value of m* /ml 

In this Letter, we will develop a non-perturbative 
statistical mechanical method to understand the novel 
strongly correlating physics. The comprehensive ther- 
modynamical quantities can be fixed by the underlying 



partition function. For instance, the entropy incorpo- 
rating the dynamical and quantum correlations can be 
derived from the grand thermodynamical potential. As 
we will see, the entropy reasonably characterizing the 
strong correlation information is crucial in determining 
the Landau parameters while ensuring the strict thermo- 
dynamical self-consistency 

The basic prospective is that the single particle spec- 
trum will be modified by the collective dispersive effects. 
The medium-dependent effective action allows a natural 
implementation of the collective effects as discussed in 
nuclear many-body literature [4j,|5j. One can imagine that 
the bare interaction potential is renormalized by the sur- 
rounding environment and the single particle spectrum 
will be modified. At unitarity, the medium dependence 
of the interaction strength can cancel the infrared diver- 
gences and lead to a finite physical result. Fixing the 
Landau Fermi-liquid parameters of strongly interacting 
fermions at unitarity is the motivation. 

In order to account for the unusual fluctua- 
tion/correlation effects, we work with the medium-scaling 
Hamiltonian in terms of the grand canonical ensemble 
framework! 16] 



H 
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In equation ([I]), a, [3 =] (a), J, (b) represent the hypcrfinc- 
spin projection Ising-variables while the m is the bare 
fermion mass. The bare coupling constant Uq = 4ira/m 
has been substituted by a medium dependent functional 
through a specific transformation 
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The n is the particle number density and [x* the effective 
chemical potential defined below through equation ([!]). 



The minus sign in the denominator of equation ([2]) 
incorporates the alternating frustration function of the 
surrounding environment. The frustration spirit is con- 
sistent with the general Le Chatelier's stability princi- 
ple in thermodynamics [16l Il7j. The nonlocal correlation 
physics characterized by the alternating minus sign coin- 
cides with the second law of thermodynamics. Further- 
more, the many-body correlation effects are taken into 
account as a spontaneously generated single-body poten- 
tial ex n r [n, T]. Sticking to the medium dependence of 
C/*rp, the additional compensatory term SH oc fi r [n, T]N a 
can enforce the energy-momentum conservation law. 

Two steps will be taken to give the grand thermody- 
namical potential with equation ((T|) and the medium- 
scaling interaction, equation p]). Firstly, the shifted rel- 
ative minimum n(T,fl) is obtained by fixing the inter- 
action analogous to the Hartree approximation with the 
linear bare potential. Secondly, the absolute minimum 
fi(T, fi) for the given chemical potential fj, is derived with 
the constraint Legendre transformation correspondence 
relations between \x and N or T and fffllilla]. 

The achieved equation of state/grand thermodynami- 
cal potential can be presented as the coupled parametric 
equations formalism 
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with 



2T 
P ideal = ~Ja h/%( z )> 



(5) 



being a function of the effective chemical potential. The 
employed notation a„a is defined by U*rr = 4na e cc/m. 

With the implicit collective variable fj,* introduced by 
the single particle Green function self-consistent equation 
(|4]), the total number density n = n-j- + nj, = 2n^ is 



n(T,„*)=2J f k = ^f 3/2 (z'), 



(6) 



with the momentum integral symbol J fc = J"<i 3 k/(27r) 3 
and quasi-particle Fermi-Dirac distribution functions de- 
fined as 



fk 



1 



z'- 1 ^ 



z' = eP"' 



1 



(7) 



The A = \j2ir /(mT) is the mean thermodynamical de 
Broglie wavelength. The (3 is the inverse temperature 
(with ks = h = 1 throughout the letter). The effective 
fugacity z' is analogous to fugacity z = e^ M . 

The first two terms of equation $5§ have exactly the 
same structure of the mean-field theory with fixed lo- 
cal interaction, which is linear-like. The last correction 
term oc n 3 is picked up in a thermodynamical way by re- 
laxing the medium-dependence of interaction potential; 



it describes the non-Gaussian correlation physics and of 
non-local characteristic. They are the rearrangement ef- 
fects of single particle energy spectrum shift reflected by 
the Brueckner-Bethe-Goldstone (BBG) technique with 
medium-dependent interaction potential [5J. From the 
statistical field theory viewpoint, this term guarantees 
the energy-momentum conservation law. 

Associated with U*rr[n,T], the shift strength C can be 



identified to be 




C(T,n*) -. 
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Furthermore, the entropy density s = S/V can be derived 
with the Hamiltonian-Jacobi implicit variable method 
from the coupled equations, © and Q 
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There is a dual relation between C and T>, which can 
be seen from equations © and (fTTj) . It is worthy noting 
that the correlation factor T>(T, /j,*) explicitly character- 
izes the temperature fluctuation in addition to density 
fluctuation. For the non-relativistic scenario, the factor 
T> can be furthermore simplified 






(12) 



Correspondingly, the entropy and energy density e = 
E/V = —P + fin+Ts can be also reduced to the compact 
formalisms with the explicit ex n 2 ,n 3 correction terms 
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• ' — - ideal - .3 -ideal' the Virial Theorem PV = \E of 
ideal gas is satisfied at unitarity[8], which is obvious by 
comparing equation and equation (|14|) . The vanishing 
ex n 2 , n 3 terms in the entropy density at T — ► ensures 
the third law of thermodynamics as expected. 

The thermodynamical properties can be characterized 
by the coupled parametric equations ([3]) and Q or equa- 
tions (|13p and (TT4")) . Relevant formulae can be expressed 



explicitly in terms of s = S/V etc. and according to T 
and the implicit variable \i* . 

For example, the sound speed squared is given by the 
adiabatic compressibility n$ — (dn/dP) s according to 



c 2 = (^-) = 1 
\dmnj g m 



dP 
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The specific heat per particle at constant volume is 

T 



Inserting the intermediate result, equation (122)) , into 
the expansions such as equation (fT9"| of the various Fermi 
integrals fj(z') (j = — |, — |, • ■ ■ , |) in Cy, c 2 and P etc., 
the thermodynamical quantities can be expanded accord- 
ing to T/T f . Here, T f = k 2 f /(2m) is the "Fermi" char- 
acteristic temperature. The calculation is quite lengthy 
by not complicated and the explicit expressions can be 
obtained. 

For example, the expansion for pressure is 



Cy 

N 



ds_ 

dT 



(16) 



The specific heat Cp per particle at constant pressure is 
calculated according to 



Cp 

N 



T 
n 






dT 



(17) 



The general expressions for c 2 , Cy , Cp etc. at finite tem- 
perature involve the high order susceptibilities described 
by such as f -5/2(2')- They are quite lengthy and are 
emitted for brevity. Obviously, the important role has 
been played by the entropy incorporating the dynamical 
high order and quantum fluctuation/correlation effects. 

Now we take the Sommerfeld lemma to do the low 
temperature expansion in order to give the Landau pa- 
rameters. To characterize the strong correlation physics, 
the expansion is in terms of T/n* — (lnz') _1 instead of 
directly according to (lnz) -1 . 

Let us start from the standard Fermi integral 
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where the factors not in the brackets are those for the 
ideal Fermi gas. By keeping up to the lowest order, equa- 
tion (|2"5|) gives the universal coefficient £ = | to which 
much attention has been paid. At unitarity, the univer- 
sal coefficient of the next order is ip It coincides with 
the first order coefficient of Cy/N, equation (|28|) . when 
\a\ = 00. However, the general expressions are quite dif- 
ferent from each other. 
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In the strong degenerate regime with /i* 3> T, 1(a) and 
Fermi integral f '5/2(2') can be approximated by 
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From equation (|19p . the low temperature expansions of 
the remaining Fermi integrals can be obtained by using 
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FIG. 1: Universal specific heat: The above solid curve is 
Cv/N of the unitary gas while the dashed one is C v /N for 
ideal Fermi gas. The labeled (a) curve is the Cv/C v ratio 
while the indicated (b) is for d 2 fj,/dndT. 

The sound speed squared is expanded as 



and the particle number density n = A:5/(37r ), the ef- 
fective chemical potential can be inversely expanded in 
terms of the "Fermi" momentum kf 
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The lengthy second expansion coefficient is still -^ at 
unitarity. By comparing with those of Landau Fermi- 
liquid theory[2|, the first expansion coefficient is 
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The lowest order of the low temperature expansion for 
specific heat gives the effective mass definition 
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The explicit difference between Cy/N and Cp/N is from 
the second order. At unitarity with \a\ — oo, the effec- 
tive mass ratio is m* /m = -tt, which agrees with the 
previous attempts 1.04 - 1.17 [HE1 El- The effec- 
tive fermion mass is the dynamical "higher order" effect 
because the linear term oc fc/a is absent in the weak cou- 
pling expansion [6||. The reciprocal term oc l/{kfd) does 
appear in the strong interaction regime. 

The universal property of unitary Fermi gas can be also 
explored by such as the specific heat at unitarity. The 



for Cy to C v of ideal Fermi gas 
^ at T — > 0; it must be saturated as Cy/C v — ► 1 
in the Boltzmann limit. It is not a naive monotonously 
decreasing function with the increase of the rcscaled tem- 
perature T/Tf. The numerical study indicates that this 
ratio has a knee peak at T* sw 0.217/ although Cy jN it- 
self is a monotonously increasing function (different from 
the finite interaction scenario) . This can be attributed to 
the competition between the particle number and tem- 
perature fluctuations. To characterize this unusual be- 
havior, the mixing susceptibility d 2 fi/dTdn is shown in 
FigfTJ where a similar knee peak around T* appears. 

To conclude, the medium scaling spirit is realized by a 
nonlinear transformation for the effective interaction and 
consequently for the physical chemical potential. The 
introduced auxiliary implicit variable-effective chemical 
potential characterizes the correlation physics and makes 
the exact grand partition function appear as the highly 
nonlinear coupled parametric-equations. 

The low temperature expansion in terms of (lnz') _1 
(consequently eliminating z') can be performed ana- 
lytically with the Sommerfeld method. The effective 
fermion mass and sound speed squared have been ex- 
plicitly calculated. As the important dynamical param- 
eter in Landau Fermi-liquid theory, the calibrated ef- 
fective mass is reasonably consistent with some Monte- 
Carlo attempts. The parametric-equations formulation 
can gauge the infrared singularity manifested by the low- 
energy long wavelength thermodynamics of the three- 
dimensional contact interaction fermions system. 
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